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Abstract. We investigate theoretically electron dynamics under a VUV attosecond pulse 
train which has a controlled phase delay with respect to an additional strong infrared laser 
field. Using the strong field approximation and the fact that the attosecond pulse is short 
compared to the excited electron dynamics, we arrive at a minimal analytical model for the 
kinetic energy distribution of the electron as well as the photon absorption probability as a 
function of the phase delay between the fields. We analyze the dynamics in terms of electron 
wave packet replicas created by the attosecond pulses. The absorption probability shows strong 
modulations as a function of the phase delay for VUV photons of energy comparable to the 
binding energy of the electron, while for higher photon energies the absorption probability 
does not depend on the delay, in line with the experimental observations for helium and argon, 
respectively. 
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1. Introduction 

Technological advance has made it possible to expose atoms and molecules to a combination 
of attosecond pulse trains (APT) and infrared (IR) laser pulses with an accurate control of their 
phase delay [1]. The photo electron spectrum of atoms in this combined light field has been 
studied [2], as well as above threshold ionization [3], the latter together with high harmonic 
generation in the combined field also theoretically [4], along with another quasi-analytical 
formulation [5] and a fully numerical R-matrix calculation [6]. 

While many parameter combinations are possible, a dynamically very interesting regime 
emerges when the energy of the VUV photon from the APT is comparable to the ionization 
potential but the IR pulse alone (typically 780 nm wavelength) is not intense enough to ionize 
the atom. The combined action of both fields leads to a time-dependent wave packet dynamics 
which is very sensitive to the phase delay, equivalent to the carrier-envelope phase (CEP) of 
the IR field. Consequently stroboscopic measurements of interfering electron wave packets 
which overlap in this energy regime but leave in opposite directions have been demonstrated 
to provide a sensitive tool to measure phase differences [7]. In a subsequent experiment it has 
been recently shown that the phase delay does not only influence the photo electron spectrum, 
but also modulates the total ionization and absorption probability [8]. A solution of the one- 
electron time-dependent Schrodinger equation (TDSE) with a pseudo potential for helium 
yields excellent agreement with this experiment [8]. While quite generally interference of 
wave packets must be responsible for the pronounced oscillatory behavior of the absorption 
yield as a function of the phase delay, the exact reason and systematics of these oscillations is 
difficult to identify in a fully numerical solution. 



CONTENTS 



3 



Here, we formulate a minimal analytical approach. It elucidates the mechanism behind 
the pronounced structures in the electron observables as a function of phase delay in the 
spirit of the "simple man's approach" successfully formulated for high harmonic generation 
in intense fields [9]. For a recent refinement of the simple man's approach, see [10]. 



2. Ionization with an APT in the presence of an IR field: Analytical approach 

We consider an electron bound with energy si exposed to a train of attosecond pulses in 
one dimension. The intensity is so weak that we are in the single photon regime for each 
atto pulse of width cr and central time Consequently we neglect the slowly varying 
envelope of the ATP and write for the interaction potential (in velocity gauge and rotating 
wave approximation) 

A''atto 

Vatto(0 = P exp(-m ^ exp[-(r - tnf/(2(r^)] , (1) 

n=l 

where p is the electron momentum and Q the central frequency of the attosecond VUV pulses. 
The infrared field is characterized by the vector potential 

A(t) = Ao cos(ajr + (p) . (2) 

It has linear polarization in the same direction as the atto pulses and a phase delay with 
respect to the APT. Its envelope Aq is assumed to be constant over the length of the APT. We 
use atomic units throughout the paper unless stated otherwise. 

Since the APT gives only rise to single-photon absorption events, first-order time- 
dependent perturbation theory provides an accurate description under which an initial state 
l^(ri)) evolves into the state at time tf according to 

mtd) = -i \' U{tu t)V,Ut)U{U h)mh)) dt (3) 



with the time evolution operators U for the electron under the combined influence of the IR 
field and the Coulomb potential of the ion [14]. 

The key idea behind the simple man's simplification which permits an analytical 
treatment is to consider the phases in the integral of Eq. (3) as the dominant contributions, 
Ua{t" , t') = exp(z[0(j,(?") - 0(j,(f')]), and to approximate the propagators before and after the 
photoabsorption differently, U{t, ti) = Ui and U(ti, t) = Uf. Then, one gets from Eq. (3) 

m = Yu r'exp(0„(O)^/?, (4) 

where the phase is given by 

(Pnit; p, if) = /[Offe) - OKO] - a - tnf/{2(r^) - i^t + m,{t) - Oi(?i)]. (5) 

Here, Oj is the phase accumulated before the VUV absorption, where the electron remains in 
its initial state, almost unperturbed by the laser field. Hence, we may write ^i{t) = -sit, and 
define the energy Sf = Si + Q. after the absorption of the VUV photon at time t'. For t > t' 
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we assume that the electron dominantly feels the IR field, described by the corresponding 
classical action (or, equivalently, quantum Volkov propagator) 



p,ip) = - ^ (/> + Ao cos(6[;t+v?))^/2 Jr. 



(6) 



A wavefunction expressed as an integral over time with two phases Oi f , as in Eq. (4), appears 
also in the analytical description of high harmonic generation, see e. g. Eq. (22) in [1 1]. In this 
context, often a stationary phase approximation is invoked. This is not possible for Eq. (4), 
since the attosecond pulse restricts the time integration effectively to an interval of the order 
of cr < In I CO, the period of the IR field. However, here we may expand the phases in a Taylor 
series about the atto peak at time This converts Eq. (4) for ti -oo and tf +oo into a 
Gaussian integral of the form 



U —CO 



k=0 



dt . 



(V) 



where the wavefunction depends now on the asymptotic electron momentum p and the phase 
delay (p. The solution to Eq. (7) reads 



In 



-F 



(8) 



The expressions are explicitly given in Appendix A. The photo-electron spectrum 
dP{p, (p)ldp and the total photo absorption probability P{ip) read in terms of 4'°°iP^ 



dP{p,<f)ldp = \r{p,<ff 



P{<p) 



I 



Jp|«A'"(p,(^)p 



(9) 



For simplicity, we will consider an absorption probability, normalized by the number of atto 

pulses A^atto 

Piif) 



(10) 



3. Explicit form of the replicated electron wave packet 

As we will see it is possible to factorize the replicated electron wave packet (EWP) into one 
term which depends on A'', the number of IR cycles over which the APT extends, and v, the 
number of atto peaks in each IR cycle. Experimentally, both y = 1 and v = 2 have been 
realized [15]. The total number of attosecond pulses is then A^atto = Nv. 

We start with the fundamental APT with one attosecond pulse in each IR period, therefore 
t„ = Inn/o). Any offset in time can be absorbed in the definition of the phase ip. To obtain 
(ff'^ip, (f) from Eq. (7) explicitly we have to evaluate the functions ^ in Eq. (8) at times 
as detailed in appendix Appendix A. Collecting all the phases from Eq. (8) but the prefactor 
{2nl{-F^n^)yi^, which contributes only logarithmically to the phases, we may write 

N-\ 

i^rip,<p) = a:J]e'^"'"^X(P,V), (H) 

n=0 
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Figure 1. Comb function K^ip) corresponding to = 8 (solid), N - A (dashed) and N - 2 
(dotted-dashed) IR cycles, for both one or two attosecond pulses per IR cycle. The IR intensity 
is / = 1.3 X lO'-' W/cm- and the excess energy Sf = 0.144 au. 



where ay is an overall phase which will not affect the observables in Eq. (9) and 

C = (//2 + t/p-£f)/a;, 
with the ponderomotive potential = A^/A. The wave packet takes the form 



(12) 



X(p, if) = exp(zpxy) exp <^ - 



(1 - i(T^o?x^{p + p^)) \ , (13) 



lalip) 

where only real valued parameters have been used. The two parameters 

Py=AoCOS^, = (Ao/6o»)sin^, (14) 

characterize the motion of an electron released at ? = into the IR field: It will quiver around 
the position + p^t having a drift momentum p^. The wave packet in Eq. (13) contains 

crAp) = 1/cr^ + [cro?{p + p^)x^) , (15) 

which is an effective width in energy with two contributions: The first one is 1/cr^, the 
variance in energy due to the temporal width of the Gaussian attosecond pulse. The second 
one accounts for the change of p^ gained from the IR field during the VUV photo ionization. 
This change is proportional to the electric field, or to x^. 

Next we evaluate Eq. (8) for two atto pulses per IR cycle, v = 2. Apart from pulses 
at tn = 2nn/oj we have a second sequence at tn+1/2 = 2(n+l/2)n/a) = {2n+l)n/a). A little 
thought reveals that for the first sequence <l>f(?„; p,(p) = Of(0; p,(p), whereas for the second 
one <l>f(?„+i/2; p, (f) = Of(0; -p, (p), see also appendix Appendix A. Collecting again all terms 
from Eq. (8) we can write iA^(/7, (p) in the form 



N-l 



(16) 



Obviously, the sum over n is the same geometric series as in Eq. (11). We call its absolute 
square the comb function. 



N-l 



n=0 



ilnirC 



sin^ (.ttC) 



Kn{p). 



(17) 
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Delay (tt) Delay (7t) 

Figure 2. Results for 1 atto per IR cycle, / = 1.3 x 10'^ W/cm^, Sf - 0.144 a.u. and = 4. 
The attosecond pulses have a width FWHM=370 as (cr =9.186). (a) <t>\{p, (f). The white lines 
represent ±(2ef)'^^ - p^. (b) Comb function Kiiip). (c) Photo-electron spectrum, 
(d) Normalized absorption probabiUty P^{(p). 

and show it for increasing A'^ in Fig. 1. From v = 1,2 the general structure of the final 
wavefunction for an arbitrary number of atto pulses v per IR cycle emerges: It 

factorizes in the comb amplitude which depends on the number of of IR cycles, and a 
complex wave packet containing v sub-packets which are created by atto pulses during one 
IR cycle and therefore depend on the phase difference (p of the IR pulse and the APT. 

4. The photo-electron spectrum 

The product structure of the asymptotic wavefunction iff'^ carries over to the photo-electron 
momentum distribution Eq. (9) since 

dPy/dp = \lfi7(p,ipt = KM^viP,^). (18) 

The function Kn(p) has maxima separated in energy by the IR frequency to which become 
sharper with increasing A^, as can be seen in Fig. 1. Hence, K^ip) acts like a comb in 
momentum for the photoelectron spectrum. The comb selects particular values of p occurring 
at specific phases ip from the electron momentum distribution Xy{p, (f), which builds up from 
the V atto pulse wave packets within one IR period. 
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Figure 3. Same as in Fig. 2 (N - 4), but for v - 2 atto pulses per IR period (A^atto = 8). The 
white hnes in (a) correspond to different branches in the wave packet X, see text. 



4.1. One atto pulse during an IR cycle 
For V = 1, we get from Eq. (13) 



[{p + p^f-l2-Sf\ 



X,{p,cp) = [y(p,y)|^ = exp|- -j;^^^^ 1, (19) 

The electron distribution Xi shown in Fig. 2a has two branches centered about p+ = 
+ V2ef - Py;. Each of them traces the streaking momentum p^ (Eq. (14) and white lines 
in the figure) which is imprinted when the attosecond pulse excites the electron with a phase 
delay if. The width cr^ of the branches has maxima at ^ = ^n, \n and minima at ^ = 0, n. 
The multiplication of X\ (Fig. 2a) with the comb K^{p) (Fig. 2b) gives the photo-electron 
momentum distribution shown in Fig. 2c. One clearly sees a preference of momenta and 
phase delays. The modulation in phase delays survives upon integration over p in the total 
absorption probability /'^'"'°((^) shown in Fig. 2d, with maxima at ^ = 0, n. 



4.2. Two atto pulses during an IR cycle 

While the comb function remains the same, we have now a more complicated single cycle 
momentum distribution composed of two wave packets during each IR cycle, 

X2{p, if) = le^'^^^'^xiP, <P) + e'^'^l^xi-P, ■ (20) 

The basic structure with two branches for each wave packet x is the same as for y = 1, 
resulting in a total of four branches at = + V2ef - p^ and p~ = + ^[2e'^ + p^, indicated as 
white lines (solid and dashed, respectively) in Fig. 3a. In addition, the wave packets interfere 
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(b) 




leading to a rich pattern in X2, as can be seen in Fig. 3a. However, again the comb K4(p), cf. 
Fig. 3b, selects specific momenta and phases for the photo-electron momentum distribution 
Fig. 3c, which produces a modulation in the total absorption probability (Fig. 3d) similarly as 
for V = 1 , with maxima at (/? = 0, n. Note, that for both v = 1 and v = 2 the number of maxima 
in the absorption probability Py((p) is the same. 

The effect of an increasing number A'^ of IR cycles in the comb K^ip) on the photo- 
electron spectrum is shown in Fig. 4 for the same intensity as before, but for an excess energy 
Sf = -0.028 and v = 2. The narrower lines for larger N (Fig. 4b as compared to Fig. 4a) is 
due to the sharper comb for larger A^. 

For the cases shown in Figs. 2, 3 and 4 the maxima appear at ^ = 0, tt and the absorption 
probabilities have similar shapes. This is not always the case. Rather, the position of the 
maxima and the contrast between maxima and minima depend on the particular comb Kn and 
therefore on N, Up and Sf, as well as on the details of the branches in the wave packet Xy. 
This will be discussed in the next section. 

5. Position of the maxima and minima in the absorption probability 

We have seen how the oscillations in the absorption probability arise from the interplay 
between the comb function K^ip) and the momentum distribution Xy of a multi-component 
wave packet. Both depend in a complex manner on the parameters of the laser fields and the 
groundstate energy of the atom. Hence the question arises if one can predict analytically for 
a given IR field intensity, where the maxima in Py((p) appear as a function of phase delay for 
different VUV photon energies Q of the APT. 

From the structure of the comb as discussed in Sect. 4 one can directly conclude that the 
oscillations in Py((p) will disappear for increasing Sf = £[ + Q, as illustrated in Fig. 5a. For 
large excess energy Sf the branches of X are centered about high absolute momentum values 
I + (2ef)^^^ ±P^\, where the comb is dense. Hence, the comb traces X(p, (p) homogeneously for 
all (f and the absorption probability hardly depends on (p (Fig. 5b). The physical meaning of 
this is that when the electronic wave packet triggered by an atto pulse leaves the nucleus with 
a high kinetic energy, the overlap with the EWPs released by subsequent atto pulses vanishes. 
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Figure 5. Photo-electron momentum distribution (left) and corresponding absorption 
probability (right) for v = 1, / = 1.3 x lO'^'W/cm^, = 8 and Sf - 0.28 a.u. (upper panels) 
andfif - -0.028 a.u. (lower panels). 



which diminishes the intereference among the wave packets. 

A more systematic analysis of the position of the maxima and minima in the absorption 
probability can be carried out analytically using symmetry properties of Py and eventually a 
stationary phase approximation with respect to p. The condition for extrema is dPy/d(p = 0, 
which can be written for the case of v = 1 as 

dP r 

-^ = J dpKMX[(p,ip) = 0, (21) 

where X\(p,(p) = dXi(p,(p)/d(p, which is proportional to x^. Therefore, the absorption 
probability has extrema for = 0. The symmetry of the functions under the integral reveals 
another set of maxima: the comb K^ip) is even in p, and for p^ = 0, the function X[ is odd in 
p. Hence, the integral in Eq. (21) is zero and there are also extrema for p^^ = 0. To summarize, 
Eq. (21) is fulfilled for every if = nn/2. 

To distinguish between maxima and minima we need the sign of the second derivative, 
X". To keep the derivation simple we will make now use of the stationary phase 
approximation, which is applicable since the comb K^ip) is a highly oscillatory function 
which depends only on p^. Therefore, its global stationary phase point is p = 0, which is also 
obvious from Fig. 1, and we get for Eq. (21) 
dPi 



dip 

where 



K^{Q)X\{Q,ip), (22) 



X\{Q,ip)=X,{0,ip) — 

dip _ 



(23) 
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Figure 6. Absorption probability as a function of for 1 (a) and 2 (b) atto pulses per IR 
cycle, for / = 1.3 x lO'^ W/cm^ and A^=8. Vertical white lines: position of the maxima at Sf 
in the two regions Sf < and Sf > 2Up, as predicted with the SPA: at (f-njl and 37r/2 in the 
first case and at 0, n and 2n in the second. The horizontal white lines differentiate these two 
regions. 



It can be easily shown that this function is zero for = or = as before, and has 
additional zeros at = 2ef. For the determination of maxima and minima we are only 
interested in the sign of the second derivative, which can be expressed for the three groups of 
extrema as 



sgn[^;'(0,^)U =sgn[^(2f/p-£f)] 
sgn[r;(0,^)|^^^„^^^^J =sgnkM 



< Sf < 2U„ 



(24) 
(25) 

(26) 



with g{x) = x{l + Acr^co^U^x)'&o-^UfXi{0,ip). This leaves a clear structure of maxima and 
minima for gf < and Sf > 2Up as summarized in Table 1 . For < < 2Up, all four extrema 
within 2n are minima, making the approximation not very trust worthy. Indeed, as will be 
demonstrated later, chaotic dynamics dominates this energy region rendering approximations 
problematic. In Fig. 6a one can see this structure of maxima and minima in (p additionally 
modulated in ef, where the distance between the maxima is oj. The latter is a consequence 
of the effect of the comb function K^iO), which has maxima at Cqtt = nn (with Co = C\p=q), 
which means maxima at energies ef = Up + noj. 

Moving on to the case y = 2 we investigate Piii'P) taking X2 from Eq. (20), and for the 
sake of simplicity we will use the stationary phase approximation from the beginning. This 
function has the same stationary phase point p = as before. Then we can write 

^2(0, ip) = A cos(Co;r/2)^i(0, ip) . (27) 

From Eq. (27) it is immediately clear that the v = 2 case has the same structure of maxima and 
minima with respect to the phase delay ^ (see table 1) as the v = 1 case. The only difference 




Figure 7. Quantum calculation for the absorption probability for an IR with /l=796nm, 
/ = 1.3 X lO'^W/cm^ and an APT with FWHM=5fs, and three different central energies: 
the 13* harmonic of the IR field (full line), the 15* (dashed Une) and the 17* (dashed-dotted 
line). The energy levels of the He atom are shown in the inset (full black lines), together with 
the ionization threshold of He and Ar (dashed lines). The energies for the 13*, 15* and 17* 
are also shown (red lines at the right). 



is the modulation in sf. The additional factor cos[Co7r/2] has maxima only at even multiples 
of the IR frequency gf = t/p ± 2nto, while it has zeros at odd multiples of to. Consequently, the 
distance between the peaks for v = 2 on the Sf axis is given by 2a) instead of as for v = 1, 
which can be seen in Fig. 6. 

6. Comparison with experimental results and exact quantum calculations 

The interesting dependence of the absorption probability on the phase delay (p was first 
reported experimentally and shown to be in agreement with a full numerical quantum 
calculation by Johnsson et al. [8]. In the meantime it has been confirmed by other experiments 
[13]. Instead of increasing the VUV photon energy (to vary Sf = Si + Q.), the ionization 
potential -ei was varied in the first experiment [8] by using He and Ar atoms as targets for the 
combined IR + APT field, the latter with two atto pulses per IR cycle and a central energy of 
n=23 eV. This energy is enough to ionize an electron from Ar, but not from He (see inset in 
Fig. 7). Strong oscillations in the ionization probability of He as a function of (p were found 
with maxima at<p = {n+l /2)n and no oscillations were detected for Ar. 

These results are in qualitative agreement with our analytical predictions, as shown in 
Fig. 6b. For < (as in He with f2=23eV), we expect maxima at (p = (n + l/2)n , 



energy 


(f = nn 


(f = {n+\/2)n 


Sf > 2Up 


maximum 


minimum 


Sf<0 


minimum 


maximum 



Table 1. Position of the maxima and minima in (p for different energy ranges Sf . 
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Figure 8. Contrast (Eq. (28)) of maxima and minima in the absorption probability for v — 2 
atto pulses per IR period. Straight line: analytical approach. Dashed line: TDSE calculations, 
for an APT with central energy from the 13* to 17* harmonics. The energies for the 13*, 15* 
and 17* harmonics are pointed out. The conditions are = 4 (A^atto = 8), / = 1 .3x10'^ W/cm^ 
and A =796 nm. The ionization threshold is marked by a vertical line. 



while for energies well above threshold as in Ar, we expect a flat absorption probability. To 
double check the transition from the positions of the maxima from (n+l/2)n to nn going 
from £f < to positive Sf, we have performed full numerical calculations. We use a 
three-dimensional one-electron model for the He atom^, and a classical electric field, with 
CO = 0.0572, / = 1.3 X 10"W/cm2 for the IR, and / = 10'' W/cm^ for the APT. Results 
are shown in Fig. 7 for three diff"erent APTs, centered at the harmonics 13*, 15"^ and 17*, 
respectively. One sees indeed that the contrast of the maxima gets smaller for increasing but 
still negative Sf, as it is the case going from the 13* to the 15* harmonic while for positive Sf 
(17* harmonic), there appear maxima at ^ = 0, n. 

Hence, the simple man's approach presented here provides in general a very good 
understanding and interpretation of the eff"ect a combined APT and IR field has on the 
ionization of atoms. Only for energies ef very close to the ionization threshold, |ef| « Up the 
simple man's approach is too drastic for reliable results. In this energy range, the oscillatory 
absorption probability depends sensitively on details of the electronic wave packet whose 
dynamics is highly chaotic. This is demonstrated in Fig. 8 for a sensitive observable, the 
contrast between maxima (at ^max) and minima (at (fimm) for difl'erent photon energies. 



One may question if such details of the chaotic behavior are helpful to understand the 
dynamics. Future work will show if more robust observables, such as correlation functions 
with characteristic correlation lengths and similar quantities are more suitable to characterize 
electron dynamics under the illumination of APTs and IR fields. 

i We propagated the TDSE in the atomic potential V{r) = - [1 + exp(-r/ro)] /r, with ro - 1.05 A guaranteeing 
the correct ionization potential of helium, and the combined laser field of APT and IR pulse. The envelope of 
the APT was a Gaussian of 5 fs width, the one of the IR pulse had a cos^ shape containing 20 cycles. 



A = 



P{<Pmax) - P{<Pmin) 
^(^max) + ^(^min) 



(28) 
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7. Conclusions and outlook 



We have presented a minimal analytical approach to understand the behavior of electron wave 
packets generated by an attosecond pulse train (APT) in the presence of a strong IR field in 
the framework of the simple man's approach in strong field physics, with special emphasis in 
the phase delay between the APT and IR pulses. In this approximation the photo absorption 
probability can be written as a product of a frequency comb function, resulting from the 
periodic nature of the pulses with the IR frequency, and the probability density of an electronic 
wave packet whose number of components is given by the number of attosecond pulses within 
on IR period. Only the latter depends on the phase delay between APT and IR field, while 
the former acts like an electron momentum filter. The minimal approach provides insight 
into the formation of oscillations in the photo-absorption as a function of phase delay, on 
the frequency of these oscillations and the general trend of the phase delay as a function of 
excess (or photon) energy. The minimal approach fails for excess energies comparable with 
the ponderomotive potential, where the chaotic nature of the electron dynamics renders results 
very sensitive to approximations. 



Appendix A. Taylor expansion of the phase 

The functions F^^^ = d^(pn{t)ldt^ needed in Eq. (8) are 



t=t„ 



Fr = .-f^^^i|^<iT (A.1) 



.(1) ip+A{tn)) 



F):> = r -isf, (A.2) 



1 



= i(p+A(tn))A'(tn)-—. (A.3) 
cr^ 

Their explicit values for = Inn I oj are 

Ff = Fq + ilnnC + ipx^ , (A.4) 



F\' = i lEf, (A.5) 



= -i{p + p^Wx^ - ^ , (A.6) 

where C, and p^ are defined in Eq. (12) and Eq. (14), respectively. Fq is a phase that does 
not contribute to \il/'^{p, (^)P (Eq. (9)). 

For two atto pulses per IR cycle (v = 2) we have in additional set of quantities F^^^^^j. 
They diff"er from the values of the F,f^ in Eqs. (A.4) -(A.6) only in two respects: The term 
ilnnC in F^^ has to be modified to i{ln+\)nC and in all quantities p has to be replaced by 
-p, which follows from of the Volkov propagator (6). 
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